Structured environments in solid state systems: crossover from Gaussian to 

non-Gaussian behavior 
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The variety of noise sources typical of the solid state represents the main limitation toward the re- 
alization of controllable and reliable quantum nanocircuits, as those allowing quantum computation. 
Such "structured environments" are characterized by a non-monotonous noise spectrum sometimes 
showing resonances at selected frequencies. Here we focus on a prototype structured environment 
model: a two-state impurity linearly coupled to a dissipative harmonic bath. We identify the time 
scale separating Gaussian and non-Gaussian dynamical regimes of the Spin-Boson impurity. By 
using a path-integral approach we show that a qubit interacting with such a structured bath may 
probe the variety of environmental dynamical regimes. 

PACS numbers: 03.65.Yz, 03.67.Lx, 05.40.-a 
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INTRODUCTION 

Controlled coherent dynamics of solid state devices 
has been demonstrated in recent years UM- Compared 
to other implementations, solid state qubits suffer from 
stronger broadband noise originating from sources with 
different character. The main limitation toward the re- 
alization of controllable and reliable quantum circuits al- 
lowing quantum computation is decoherence due to ma- 
terial (and device) dependent noise sources. The result- 
ing noise spectrum is non-monotonous and sometimes 
characterized by resonances. Often these features may 
be attributed to interaction with a nonlinear and non- 
Markovian environment [3| . 

In superconducting nanocircuits a particularly detri- 
mental role is played by fluctuating impurities located 
in the insulating material surrounding the qubit, which 
are responsible for charge noise and flux noise 0, 01 ■ 
Background charges are known to be responsible for 
low- frequency 1// noise moreover experiments with 
Josepshon devices suggested that spurious two-level sys- 
tems may also affect high-frequency noise 0. Connec- 
tions between low and high-frequency noise features have 
been suggested in the recent experiment Ref. Q. Differ- 
ent microscopic mechanisms [9!] and effective models 
have been recently proposed to explain the observed spec- 
tral features. 

Predicting decoherence originating from such a struc- 
tured environment often responsible for non-Gaussian 
noise is a non trivial task, which has attracted a lot of at- 
tention in the past years d, [ll|, [l^, l3, 14 1 . A well estab- 
lished scheme consists in studying the reduced dynamics 
of an extended system composed of the qubit and of the 
environmental degrees of freedom responsible for non- 
Gaussian behavior. This strategy, combined with an ap- 



propriate classification of the noise sources (i.e. adiabatic 
or quantum noise), each treated via appropriate approx- 
imate tools, provides a general scheme to deal with the 



variety of noise sources typical of the solid state |14l . 
In this paper we focus on a prototype impurity model, a 
two-state impurity linearly coupled to a dissipative har- 
monic bath. Such Spin-Boson models have been thor- 
oughly investigated with a variety of methods since the 



'80s 3) 12 1- Thus the impurity dynamics is well known 
in a wide region of parameters space. This allows the 
identification of the impurity characteristic time scales 
and, therefore, of the conditions where deviations from 
Gaussian and/or weak coupling regimes are expected. 
Having this information at our disposal we will apply 
standard techniques developed for quantum dissipative 
systems to find the qubit dynamics in the presence of 
this structured bath. The analysis will provide evidence 
for the appearance of non-Gaussian effects. In partic- 
ular, their onset will be shown to be related to clearly 
identifiable effects in the qubit behavior. 

In Section 2 we introduce the qubit-impurity model 
and identify the relevant impurity dynamical quantities. 
In Section 3 we review the equilibrium correlation func- 
tion for a Spin-Boson impurity and identify its correla- 
tion time. In Section 4 we study the qubit dynamics in 
the presence of the Spin-Boson impurity within a path- 
integral approach and discuss the main features of the 
crossover from weak to strong coupling. 



MODEL AND RELEVANT DYNAMICAL 
QUANTITIES 

To be specific we shall refer to superconducting qubits 
based on the Cooper-pair box P, [l3|. Under proper 
conditions the device behaves as a two-state system de- 
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scribed in terms of Pauli matrices by (h = 1) 



i^qubit — ^ z (-) 



(1) 



The charging energy Ec gives the additional cost of 
adding an extra Cooper pair to the superconducting is- 
land and the possibility of coherent transfer of pairs 
through the junction is given by the Josephson term 
Ej (Ta;/2. The charge on the superconducting island may 
fluctuate because of interaction with uncontrolled impu- 
rities. Here we model a single impurity with a Spin-Boson 
model, the overall Hamiltonian being given by 



n = 

Ti-SB 



e 



A 



1 



Xt,+He- 



(2) 

(3) 




FIG. 1: r2(r) of Eqs.dTTI and (fT5)l as a function of tempera- 
ture. Inset: non-monotonous behavior for T -C T* = 31 Ar. 
Parameters are K = 0.01, ojc/Ar = 31, and we fixed fcs = 1- 



The two-level-system impurity (f) is coupled to a har- 
monic bath, described by TCe = X]q ■^a'^a'^aj '^i^ 
collective coordinate X. Its effect on the impurity de- 
pends only on the spectral density G{uj) or equivalently 
on the power spectrum 



S{co) = 



dt i {X{t)X{Q) + X{{))X{t)) e"^* 



ttG{\u}\) coth 



(4) 



where (...) denotes the thermal average with respect to 
H_E, and /3 = l/ksT. We consider the standard case 
when the coupling operator is a collective displacement 



X = J2a ^a{ao 



with ohmic spectral density 



G(w) = J2xld{uj-0Ja.) ^ 2K\cj\ 



(5) 



where lUc represents the high frequency cut-off of the har- 
monic modes. 

A first step in understanding the effects of damping is 
to view the impurity f and the ohmic bath as an envi- 
ronment for the qubit a. This environment is in general 
non-Gaussian and non-Markovian. A Gaussian approxi- 
mation of this structured bath amounts to replace it with 
an effective harmonic model directly coupled to cr^ and 
with power spectrum 5'i-(w) 

SAu;) = l pdt{{TS)rM+rMrS)) " (T.)L)e^"*, 

(6) 

the Fourier transform of equilibrium symmetrized auto- 
correlation function of the impurity observable which di- 
rectly couples to the qubit. Here the thermal average is 
performed with respect to Hsb and {tz)oo is the thermal 
equilibrium value for Tz- Under this approximation and 
using a master equation approach, the relaxation and 
dephasing rates for the qubit a in lowest order in the 



coupling V read |17l , [2 
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V^SriE) 



T* 



^V^Sr{E) fEc\^V^Sr{0) 



E 



(7) 



(8) 



E^j is the qubit splitting. The va- 



where E — \/E'^ -r j 
lidity of this standard approach is limited to couplings 
V ^ I/tc [i^l, where Tc is the range of the correlation 
function in Eq.® (to be defined in the next Section). 
Clearly if the impurity r has a slow dynamics vtc S> 1 , 
this picture does not apply and we have to resort to other 
methods. However previous studies on a similar model 
have shown that the Gaussian approximation may give 
good results even for vt^ > 1 but for shorter and shorter 
times, as long as Tj, increases 0]. 

^From a different perspective the failure of the Gaus- 
sian approximation can be understood by viewing the 
qubit (T as a measuring device ^2)\ for the mesoscopic 
system described by the Spin-Boson model involving r. 
A rather rough measurement protocol (short times, av- 
eraging of results) makes the dynamics of a essentially 
sensitive only to Sr{i-^), whereas if the Spin-Boson has a 
slow dynamics the spin a is able to detect also details of 
the dynamics of t which go beyond Sr{uj), and have to 
be described with more careful methods. 

In the following we will treat the impurity t on the 
same footing as the qubit a, we will apply standard meth- 
ods to trace out the bosonic degrees of freedom without 
any approximation on the qubit-impurity coupling. 
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FIG. 2: Quality factor of the damped oscillations Q{T) — 
Q{T)/'y{T) from Eqs.([TT]), ([12} and HS]), dHJ. Parameters 
are fixed as in Fig[T] 




T/Ar 

FIG. 3: Dephasing rate y(T) from Eqs.([l2} and ^ (full 
Une), for temperatures T < T*{K) and relaxation rates 
71/2(3") from Eq.dTTll for T > T*{K). The rate 71 (dotted) 
increases with temperature, 72 (dashed) shows the Kondo be- 
havior. The value of T*{K)/Ar Ri 31.8 is indicated by the 
dashed gray line. Parameters are fixed as in Fig[T] 



IMPURITY DYNAMICS AND CORRELATION 
TIME 

In this Section we will identify the characteristic time 
scale of the dynamics of the equilibrium fluctuations of 
the Spin-Boson impurity described by Eq.Q. The un- 
coupled (i.e. for V = 0) impurity dynamics strictly de- 
pends on the damping strenght K and on the tempera- 
ture |17[ . We consider the small damping K <^ 1 regime 
where series of crossovers from under-damped to over- 
damped oscillations and to relaxation dynamics with in- 
creasing temperature are observed. This analysis is also 
relevant to understand the effect of ensemble of impu- 
rities when a wide distribution of parameters has to be 
taken into account as in Ref. In fact at any fixed 

temperature each impurity may display a specific dynam- 
ical behavior. As a consequence, various sets of impuri- 
ties may effect the qubit dynamics in qualitatively differ- 
ent ways. 

Our goal is to establish to which extent the qubit may 
probe the variety of regimes of its environment dynam- 
ics. The important scale allowing the identification of 
Gaussian/non-Gaussian dynamical regimes is extracted 
from the equilibrium auto-correlation function of the ob- 
servable Tz which directly couples to the qubit. We re- 
mark that the thermal initial state of the Spin-Boson sys- 
tem, which is implied by the equilibrium correlation func- 
tion, may originate peculiar time-dependencies. Quali- 
tatively different behaviors may in fact be displayed by 
the non- equilibrium correlation function, which is evalu- 
ated for a factorized initial state of the Spin-Boson sys- 
tem [illiii. Evaluation of equilibrium correlation func- 
tions for the Spin-Boson model is a non trivial task. How- 
ever it has been shown 17| that in the unbiased case 
e — 0, and for small damping K ^ 1, the Tz auto- 
correlation function does not depend on the initial cor- 
related or factorized state. Here we focus on this case 
and recall the characteristics of St{uj) more relevant for 
our analysis, the interested reader may find details of the 



derivation in Refs. 17. loj. 

The crossover from under-damped to over-damped be- 
havior with increasing temperature takes place at 



T*{K) 



TiKkf 



(9) 



where A,. = A(A/wc)^^*-^~^'' is the renormalized tun- 
neling amplitude in the Spin-Boson model. 

The explicit form for Sr{'^) depends on 
temperature 22] 
For T < T*(k) it reads 

7 (fi - w) tan 7 -t- (f2 -I- w) tan 
ST{i^) = — aT^tT — ^ + — , , — (10) 



{lj - ny + 72 



{lo -f + 72 



where i&ncj) = 7/^^. The effective frequency SI and the 
relaxation rate are 



= 1 



l{T) = 



Re^fi^!_) _ In f^^) - ^(1) 
. 2TikBT) \2TTkBTJ ' 

(11) 
(12) 



for temperatures fcsT < A^ and 

f2(r) 



7^ s 2-2K 
l-(^ 



7(r) - -KKkeT = - 



(13) 
(14) 



for larger temperatures A,. <kBT< kBT*{K). 

Here 2r is the white noise level S{u)) « 2r for fre- 
quencies w ^ 27r//3 in (III). Oscillators with frequencies 
27: / P < uj <^ ujc renormalize the tunneling amplitude to 



At = Ar(27rfcsT/Ar) 



K 



(15) 



In Fig[T]we show the overall behavior of ^{T). At very 
small temperatures fc^T <C A^ (see inset in FigH]) (l{T) 
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increases oc T^. The two forms ^ and ^ smoothly 
match on each other at fcsT = A^. A maximum is 
reached at T = K'^/'^i^-K)^* > A^/fes, above this tem- 
perature ri(T') decreases monotonously and approaches 
zero at T*. Coherent oscillations are dephased on a 
scale 1/7 (T) which decreases monotonously starting from 
1/7(0) = 2/ {irKAr). The resulting quality factor of 
the damped oscillations Q{T) — Cl{T)/^{T) is shown in 
FigH 

For temperatures higher than T* the dynamics is in- 
coherent and Sr{^) has a different form 



Sr{uj) = 2- 



72 



71 



7l/2 



72 



71 w2 -t- 72 



+ 2- 



71 



72 



71 - 72 -I- 7J 



(16) 



r\2 



(i) 



A|. 



(17) 



Note that for T ^ T* (K) one of the two rates increases 
with temperature, 71 — > T, whereas the other shows the 
characteristic Kondo behavior, decreasing with tempera- 
ture 



1-2K 



(18) 



These features are illustrated in Fig (21 

The typical scale of the equilibrium fluctuactions of 
the Spin-Boson environment described by Sr(j-o) defines 
the correlation time Tc- Usually in the literature one 
is faced with environment models characterized by dy- 
namic fluctuations which tend rapidly to zero with time. 
In these cases represents the order of magnitude of 
the width of the environment fluctuactions [20|. In the 
present case however looking at the Spin-Boson impurity 
as an environment characterized by Sr{io) the identifica- 
tion of Tc is not immediate due to the different forms of 
the equilibrium fluctuations given by Eqs.(fTO| and (flG]) . 
For high temperatures S'r(w) is approximately a single 
Lorentzian centered at a; = and width 72, leading to 
the identification Tc ~ 1/72(7), see Fig ID In the opposite 
limit of very low temperatures S'r(w) has a double peak 
structure representing a bath responsible for oscillating 
fluctuations very weakly damped. In this unusual envi- 
ronment regime, the typical scale of the impurity fluc- 
tuations is represented by l/(l{T) which plays the role 
of Tc- For intermediate temperatures in general two al- 
most superimposed Lorentzians contribute to Sriuj) and 
Tc may be approximately identified from the condition 
Sril/Tc) = 5;"'^^/2. The resuffing Tc{T), illustrated 
in FigO interpolates between the asymptotic behaviors 
at low and high temperatures. The slight reduction of 
Tc{T) at intermediate temperatures is a consequence of 
the crossover from under-damped to incoherent dynam- 
ics, as shown in FigH) 

Once the impurity correlation time is identified, the 
condition vTc(T) = 1 separates the weak-coupling regime 
of the qubit dynamics, from the strong coupling regime 
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FIG. 4: Equilibrium correlation function 5't(w) for increas- 
ing values of the temperature. Top: from under-damped to 
over-damped regime and T/Ar — 3 (light gray), T/Ar — 15 
(black). Dashed lines indicate the width 1/tc at half height. 
Bottom: from over-damped to incoherent regime T/Ar — 20 
(light gray), T/Ar = 30 (dashed), T/A,. = 40 (black). Pa- 
rameters are fixed as in Fig[l] 



where non-Gaussian behavior shows up. In the first case, 
when vTc(T) <C 1, the standard master equation predicts 
exponential decay with the Golden Rule rate ^ given in 
Eq.®. /^From the above analysis we expect the master 
equation result to be valid in the following regimes: For 
temperatures T < T* if v/d{T) « v/Ar < 1, for larger 
temperatures T > T* if w/72(T) « wF/Af, < 1, this 
condition can be cast in the following form 



a: « (2r) 



(19) 



Therefore, for small values of v/Ar a crossover from 
weak to strong coupling is expected with increasing tem- 
perature. For the ensuing discussion here we report 
the expected value of the pure dephasing rate I/T2 = 
{Ec/Efv^Sr[0)/2 with Eqs.dlOl) and ^ 



T <T* (20) 
T > T* (21) 



1 _ ^Ec\'^ 4v'^j{T) 



To* 



(at) 



r 



The two forms match on each other at T*, and it is easy 
to show that Eq.imi) approximates I/T2* also for T < T*. 
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FIG. 5: Correlation time Tc{T): for sufficiently small temper- 
atures Tc ~ 1/Ar. For T 3> T*{K) the asymptotic behavior 
~ 1/72 (r) is indicated (dashed). The interpolating form for 
intermediate temperatures has been obtained from the condi- 
tion St{1/tc) = Sr'"'/2. Parameters are fixed as in FigH] 



QUBIT DYNAMICS: PATH-INTEGRAL 
APPROACH 



The discussion of the previous Section has evidenced the 
existence of a large parameter regime where the Gaussian 
approximation of the Ohmic Spin-Boson model does not 
apply. In this Section we study the qubit dynamics via 
a path-integral approach which includes as a special case 
the regime of Gaussian behavior of the impurity dynam- 
ics. We find exact expressions which we discuss for finite 
temperatures, specifically for fcsT > y/Sf+'iP- 

We focus on the so called pure-dephasing regime, Ej ~ 
0, which represents the point of maximum noise sensitiv- 
ity of the qubit. Thus the more interesting in the per- 
spective of using the qubit as a "noise" analyzer. In the 
pure dephasing regime the charge on the qubit island is 
a constant of motion since [H, (Jz] — 0, dephasing be- 
ing described by the decay of {o^jy) or equivalently of 
the coherences (cr±). A simple analysis shows that the 
coherences are related to correlation functions involving 
the Spin-Boson variables, specifically we found [3, [3] 



('^-(0)) 



JEct 



JEct 



Trssjc 

C-+{t) 



-iHsB-t 



Pt{0) (g) Wf3 e 



(22) 



where we have chosen a factorized initial density ma- 
trix for the qubit-impurity, p(0) — PaiO) (8) Pr(0), with 
the impurity r initialized in the mixed state /3t(0) = 
^ I + ^Sp{0) Tz and the bosonic bath in its thermal 
equilibrium state Wfj. The two conditional impurity 
Hamiltonians 'Hsb± depend on the qubit state and read 

'HsB± = T^SB ± f Tz- 





0.0. 




-0.5 








—1 .U - 


CD 




DC 


-1.5- 




-2.ol 







1.0 



0.0 




-1.0 



FIG. 6: Top: real parts of the exact solutions Ai of the 
pole equation D{X) = as a function of temperature for 
ksT > y/v'-' + A'i « 0.87Ar. Bottom: corresponding imagi- 
nary parts. Two Ai are complex conjugate below T- « 6.1 
and above T+ ~ 7.2 Ar. All Ai are real at intermediate tem- 
peratures T- < T < r+. The dominant pole is real (black 
dashed) until T+ where the character of the dominant solu- 
tion changes. For T > T+ the dominant poles are complex 
conjugate. Parameters are v/Ar = 0.25, u>c/Ar = 30 and 
K = 0.05. 



form of the correlator C ^-{t) reads 

1 



C-+iX) = [ A + A'l ( A) ~ tvSpiO) ] (23) 

i:>(A) ^ + XKi{X) +ivK2{X) . (24) 



An exact formal series expression in A for the ker- 
nels Ari(A), Ar2(A) has been derived in Ref.[2^. In the 
Markovian regime for the harmonic bath, i.e. for K <^ 1 
and temperatures -^Z + < fcsT ^ uJc, all contri- 
butions to Ari(A) and Ar2(A) of order higher than 
cancel out exactly. The lowest order contributions [2^ 
do not depend on the coupling v and coincide with the 
kernels entering the dynamics of r in the uncoupled case 
(v = 0) fr?] which read 



/Ci(A) = A| 



A + r 



e2 + (A + r)2 



/C2(A) 



-ttATA^ 



(25) 
(26) 



^ e2 + (A + r)2 ■ 

In Ref. 23] it has been shown that the Laplace trans- Inserting Eqs. ipS)) - (pS)) in (|23p and (I24p C |-(A) is read 
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FIG. 7: Top: C_+(t) for v/K = 0.25 and /tsT/Ar = 4 
(black), ksT/Ar = 15 (gray). Bottom: corresponding 
Fourier transform. Here K = 0.05, ujc//^r = 30, 5p{Q) = 0. 



ily found as 



C-+(A) 



[(A+r)^+e^] [A-wjp(0)]+A^(A+r) 
Z5(A) 



D{\) = (A2 + z;2) [(A + r)2 + + A2,A(A + r) (27) 
—iTiKveA'^ . 



The scales entering the time evolution of C ^{t) are 

found from the solution of the pole equation D{X) — 0, 
which have been reported in Ref . (isj . Here we specify 
to the unbiased case e = 0, the goal being to elucidate 
the correspondence with the expected Gaussian/non- 
Gaussian dynamical regimes as deduced from the equi- 
librium correlation function S'r(ti') discussed in Section . 



Pure dephasing due to a unbiased impurity 

When e = the pole condition ^(A) = with Eq. fTT]) 
reduces to a cubic equation which has either one real and 
two complex conjugate solutions, or three real solutions. 
We denote the three roots as Ao = — Aq € TZe and A1/2 — 
—A ± iSE, where 5E is either real or purely imaginary. 

The expression of C \.{t) in terms of the A^ is obtained 

by inverting the Laplace transform (|27p and reads 



FIG. 8: Top: real parts of the exact solutions Ai of the 
pole equation -D(A) = as a function of temperature for 
ksT > \/v'^ + Ar ~ 1.2Ar. Bottom: corresponding imagi- 
nary parts. Two Ai are complex conjugate and one is real 
for any temperature. The dominant poles are always com- 
plex conjugate. In this regime no crossing takes place among 
the 7?.e[Ai] and the dominant root is non-monotonic with T . 
Parameters are n/A^ = 0.8, K — 0.05, uJc/Ar = 30. 



+B sin {6Et) e-^* 
-2AoA + A| - i2Sp{0) vA 



B 



(A- Ao)2 + (5£;2 
[A(l -A)+ AaA - i6p{0)v] 
SE 



"(28) 
(29) 
(30) 



It is possible to show that the character of the roots de- 
pends on v/At and on the temperature. In particular 
for V < At/2^/2 we can identify two temperatures 



± 



1 r /w\2 5A2 
VA8 - 24i;2A6 + 192u4A4 - 512i)6A2 



1/2 
(31) 



32^2 

such that for T < T- and T > one solution is 
real and two are complex conjugate, whereas for inter- 
mediate temperatures the three solutions are real. For 
V > Ar/2-\/2 there is always one real and two complex 
conjugate solutions. Analytic expressions for the roots 
are quite cumbersome, approximate forms have been re- 
ported in [2^. Here we discuss the physically relevant 
regimes where crossover are expected in the behavior of 
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FIG. 9: Top: C^ + (t) for v/Ar = 0.8 and fcsT/Ar = 4 
(black), ksT/Ar = 15 (gray). Bottom: corresponding 
Fourier transform. Here K = 0.05, ljc/A,. = 30, Sp{0) = 0. 



C i^{t). To this end we focus on the dominant Ai, i.e. the 

smallest in absolute value. The analysis is conveniently 
performed distinguishing regimes where v/At < 1 or 
v/At > 1 

Case v/At ^ 1 In this regime the characteristics of 
the dominant pole change qualitatively with increasing 
temperatures (with lower bound fc^T > ^yv^+A^). For 
temperatures low enough to fulfill the condition 



At 



1 



4ttKv 



kBT+ 

Ar 



(32) 



the dominant scale is real and reads Aq ~ (^)^r. In 
dSni) v4 « 1 - (r/2AT)^ and B w -iSp{0)v/AT therefore 



C+-{t) w exp [-Aot] 



(33) 



With increasing temperature, above 7+ the dominant 
scales are complex conjugate, X1/2 — ± iSE w 
— Ay/2r± iv. Figure [6] illustrates the crossover between 
the two regimes for v/Ar = 0.25. Since A <^ v/T and 
B w ~iSp{0), we get 



C+_ (t) w cos {vt) e 



-At 



i5p{0) sin {vt) e 



-At 



(34) 



As expected, for small enough temperature the Gaussian 
approximation for the structured bath applies and a 
single scale dominates the qubit dynamics. It is easily 
seen that the condition ([32|) corresponds to (fT9|l which 
was derived from the weak coupling criterion vTc ^ 1. 
Moreover Aq = l/T^ = 5^(0)72 as given in Eq. (PT|l . 



For higher temperature non Gaussian effects show up. 
The qualitative change is characterized by damped os- 
cillations at frequency f« w in C-| (t) as shown in FiglTl 

and beatings at frequencies Ec ± v in the coherences 
Eq. ip^ . Note that the coupling strength v only enters 
the induced frequency shift, whereas the decay rate 
shows the Kondo behavior A oc T^^-^^ cfr Eq. pS]) . 



Case v/At ^ 1 For larger values of v/At the system 
stays in the regime where the dominant scales are com- 
plex conjugate and show a nontrivial temperature depen- 
dence. At temperatures T > yV + A^, the two poles 
read A1/2 ~ —{At/v)'^T/2 ± iv, with increasing T in- 
stead Ai/2 ~ — A|,/(2r)±w. The qubit dynamics follows 
Eq. (p4|l . The dominant rate A is non-monotonous first in- 
creasing with temperature and than decreases oc x^^~^. 
For large enough temperatures the decay of the coher- 
ences does not depend on v, as expected. 

This qualitative behavior is already present for inter- 
mediate values of v/A, as shown in Fig[8l The poles Xi 
never cross, therefore there is no change in the character 
of the qubit dynamics. In the specific case considered 
at small/intermediate temperatures the real parts of the 
three poles are of the same order. This is reflected in the 

Fourier transform of (t) where three Lorentzians can 

be identified, one centered at u; = 0, the others at « ±v, 
see FiglH 



DISCUSSION 

The presented analysis has shown that a damped im- 
purity may behave as an effective short-time correlated 
or as a non-Gaussian environment depending both on its 
coupling with the qubit and on temperature. The rele- 
vant scale separating the two regimes is given by Tc- In 
Section the correlation time of the unbiased Spin-Boson 
model has been found for small damping K <^ 1 at any 
temperature. We have shown that the qubit may act as a 
detector of non-Gaussian dynamical behavior, the most 
evident effect being the occurrence of beatings which are 
expected with increasing temperature. We remark that 
the results we have illustrated have been derived within 
the NIBA which limits temperatures to values larger than 
^/^P + A^. A interesting issue is to analyse the small 
temperature regime where we expect that crossover ef- 
fects may take place also for large values of v/At- A 
detailed analysis of the low temperature regime can be 
performed within the systematic weak damping approx- 
imation [l3| and will be reported elsewhere 24 1 . 
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